Dislocation Free Island Formation in Heteroepitaxial Growth: An Equilibrium Study 
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We investigate the equilibrium properties of strained heteroepitaxial systems, incorporating the 
formation and the growth of a wetting film, dislocation free island formation, and ripening. The 
derived phase diagram provides a detailed characterization of the possible growth modes in terms of 
the island density, equilibrium island size, and wetting layer thickness. Comparing our predictions 
with experimental results we discuss the growth conditions that can lead to stable islands as well as 
ripening. 
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Heteroepitaxial growth of highly strained structures 
has gained interest lately as it offers the possibility to 
fabricate nanoscale islands with very narrow size distri- 
bution [fij. Typically, H monolayers of atoms are de- 
posited on a substrate, where the substrate and the de- 
posited film have different equilibrium lattice constants. 
For small coverage the experiments document the pseu- 
domorphical formation of a wetting film, but after the 
film reaches a certain critical thickness, H Cl dislocation 
free islands form on the substrate. Thanks to their small 
and uniform size, these islands, coined self-assembling 
quantum dots (SAQD), are candidates for three dimen- 
sional electron confinement 

The controlled production of SAQDs for both opti- 
cal and electronic applications requires a good descrip- 
tion of the basic mechanisms determining the size and 
the distribution of the islands. However, such an under- 
standing is hampered by the coexistence of equilibrium 
and nonequilibrium effects: while the experimentally well 
documented existence of a flux independent critical wet- 
ting film thickness H c , is well described by equilib- 
rium theories of heteroepitaxial growth ||] , the observed 
flux and temperature dependence of the island sizes (!]] 
provide direct evidence of nonequilibrium effects con- 
tributing to the island formation process Q]. A detailed 
theory of SAQD formation, incorporating both nonequi- 
librium and equilibrium effects is beyond reach at this 
point. However, since for reversible systems nonequilib- 
rium effects represent the path of the system towards an 
equilibrium state, an adequate theory of SAQD forma- 
tion should first provide a detailed description of the equi- 
librium states supported by the dislocation free strained 
system. An important step in this direction was taken by 
Shchukin et al. |^] , who found that depending on the ma- 
terial constants and the misfit one can obtain either sta- 
ble islands, or ripening takes place in the system. How- 
ever, by neglecting the existence of the wetting layer, 
their study could not predict the actual growth mode, 
nor could provide the island density, island size, and the 
wetting layer thickness as a function of the deposited ma- 
terial, quantities that can be measured experimentally 
with great accuracy jj]. 



In this paper we investigate the equilibrium proper- 
ties of strained heteroepitaxial systems, incorporating the 
growth of the wetting film, dislocation free island forma- 
tion, and ripening. Our results can be summarized in a 
phase diagram, that not only predicts the main growth 
modes, but also provides a detailed characterization of 
the possible phases in terms of the island density, equi- 
librium island size and wetting layer thickness. We find 
that the stability of the islands depends very sensitively 
on the coverage, i.e. the misfit strain and the coverage 
have to exceed a critical value for stable islands to exist, 
and that for any misfit there is a second critical coverage 
beyond which ripening occurs. 

Model and free energy — We consider that H mono- 
layers of atom A with lattice constant dA are deposited 
on top of the substrate B with lattice constant ds, and 
are allowed to equilibrate. Due to the lattice mismatch, 
e = (dA — c?s)/rfB, in equilibrium one expects that a cer- 
tain fraction of the atoms A forms a wetting film of n\ 
monolayers and the rest of the material (H — n\ mono- 
layers) is distributed in 3D islands. We consider that the 
3D islands have a pyramidal shape with a fixed aspect 
ratio, corresponding to a single minima in the Wulff's 
plot. Neglecting evaporation, the deposited material rep- 
resents a conserved system in equilibrium with a thermal 
reservoir, thus the relevant thermodynamical potential 
density is the free-energy per atom, / = u — Ts, where u 
is the internal energy density, T is the temperature and 
s is the entropy density of the system. However one can 
show that the entropic contribution to / is negligible, 
thus / « u, where 

u{H,ni,n 2 ,e) = E m i(rii) + n 2 E is i + (H - ri\ - n 2 )E rip . (1) 

The first term provides the contributions of the ri\ 
strained ovcrlayers, being an integral over the binding 
and the elastic energy densities. The energy density of 
a uniformly strained layer is given by G = Ce 2 — $>aa, 
where — $>aa is the energy of an AA bond and C is a 
material constant, being a function of the Young's mod- 
ulus and the Poisson's ratio [||. At the wetting layer- 
substrate interface atoms have AB bonds with the sub- 
strate with a binding energy —&ab, such that A = 
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*f>AA — $ab < (wetting condition). However, due to 
the short range intermolecular interactions the binding 
energies of A atoms close to (but not at) the substrate is 
also modified [0: as we move away from the substrate, 
the binding energy density increases from —<&ab (in the 
first monolayer) to its asymptotic value — &aa- These 
intermolecular forces are responsible for the critical layer 
thickness larger than one monolayer in heteroepitaxy || . 
To include this effect we calculate the total energy stored 
in the wetting layer as 

E m i{n 1 )= / dn{G + A[(e(l-n) + e(n-l)e-^/ a ] 
Jo 

where Q(x) = if x < and 6(x) = 1 if x > 0. The 
a = limit corresponds to the absence of the short range 
forces. While (|J) provides a reasonable fit to the result 
of Ref. H , the particular form of (^) does not modify the 
qualitative behavior of the free energy provided that the 
binding energy is strictly monotonous and bounded as a 
function of n. 

The second term in eq. (|l|) describes the free energy 
per atom of the pyramidal islands and the island-island 
interaction Q 

E isl = gCe 2 - <P A a + E Q (~ lne 1 /^ + SL + ^) , ( 3 ) 

V x A x x i / z / 

where x = L/Lq is the reduced island size, Lq being 
a material dependent characteristic length [^). Depar- 
tures from planar geometries can lead to the relaxation 
of the strain energy. Thus the strain energy density of 
the islands (first term in eq. (||)) is lower than that of 
the compressed wetting layer, this reduction being ex- 
pressed by the form factor g (0 < g < 1) g. The second 
term stands for the binding energy. The elastic energy 
of an edge of length L is proportional to — LlnL [Q, 
thus the energy density is ~ — InL/L 2 , accounting for 
the first of the three terms in the parenthesis. The in- 
teraction of the homoepitaxial and heteroepitaxial stress 
fields leads to a cross term — e/L ~ —e/x. Furthermore, 
the facet energy is proportional to the area of the facet, 
L 2 , giving the energy density as ~ 1/L ~ 1/x. The 
cross term and the facet energies are combined in the 
second term in the parenthesis of eq. (3), a/x, with 
a = p(7 — e), where p and 7 are material constants de- 
scribing the coupling between the homoepitaxial and het- 
eroepitaxial stress fields (also function of the island ge- 
ometry) and the extra surface energy introduced by the 
islands, respectively. Finally, since the stress fields of the 
individual islands overlap, there is island-island interac- 
tion, described by the last term in the parenthesis of eq. 
(|), where (3{n) = be 2 n 3 / 2 @. This can be expressed 
in terms of the average island spacing d = 1/ y 'pi S i and 
the reduced island size x, giving the interaction term as 
~ (x/g?) 3 , corresponding to the dipole-dipole interaction 
between the islands ||]. The energy terms appearing in 
eq. (3) are scaled by the characteristic energy Eq set 
by the edge energy of an island of size Lq. We also scale 



C, <f>AAj and <&ab by Eq, thus the results are independent 
of the numerical value of Eq. 

The total elastic energy density of the ripened islands 
can be obtained from (|3|) by taking the limit x — > 00, 
providing E r i p = gCe 2 — &aa, which is multiplied by 
the total number of atoms stored in the ripened islands, 
(H-m- n 2 ). 

Phase diagram — Eqs. (|l])-(||) define the free energy of 
the wetting film and 3D pyramidal islands, whose min- 
ima determines the equilibrium properties of the system. 
Consequently, we have to minimize / in respect to m, 112, 
an<£l2^. The growth modes (phases) provided by the min- 
imization process, as a function of the two most relevant 
experimental parameters, the amount of the deposited 
material H and the misfit e, are summarized in the phase 
diagram shown in Fig. 1. In the following we discuss the 
properties of the phases predicted by our analysis. 

FM Phase: If H < H Cl (e), the deposited material con- 
tributes to the pseudomorphic growth of the wetting film 
and islands are absent, reminiscent of the Frank van der 
Merve (FM) growth mode. The free energy has its min- 
ima at n 2 — and m = H, indicating that the thickness 
of the wetting layer coincides with the nominal thickness 
of the deposited material, H (see Fig. 2a). The growth 
of the wetting layer continues until H reaches a critical 
value, H Cl (e), the phase boundary between the FM and 
the i?i or SKi phase. 

i?i Phase: Above H Cl (e), for < e < e±, the free en- 
ergy develops a new minima at n 2 — and < n\ < H . 
Consequently, after the formation of a wetting layer of 
m = H Cl (e) monolayers, the excess material (H — n\) 
contributes to the formation of ripened islands. The free 
energy decreases monotonically for large x, thus there is a 
tendency to accumulate all available material (H — H c (e)) 
in as large islands as possible. These ripened islands, be- 
ing infinitely large, have zero density. 

SKi Phase: Above H Cl (e), for t\ < e < e 2 , the free en- 
ergy develops a new minima at nonzero ri\ and n 2 , such 
that n\ + n 2 = H, i.e. the deposited material (H ) is dis- 
tributed between n\ layers, forming the wetting film, and 
finite islands, whose total mass is n 2 (see Fig. 2a), similar 
to the Stranski-Krastanow (SK) growth mode. At H Cl (e) 
the equilibrium island size jumps from zero (in the FM 
phase) to some finite xq(H,6) value. Naturally, within 
the SKi phase, ni, n 2 , xq, and the island density p are 
continuous functions of H and e. As Fig. 2a indicates, 
with increasing H, p increases from zero at H c to a finite 
value. Due to the island-island interaction, the wetting 
layer continues to grow sub-linearly. 

R 2 Phase: In this phase, the free energy surface has 
a minima at < n\ < H, < n 2 < H, such that 
H — ni — n 2 > 0, indicating that the deposited mate- 
rial is distributed between a wetting film, finite islands, 
and ripened islands (H — n\ — n 2 ). The finite islands 
formed in the SK^ phase is preserved (see Fig. 2a and 
b), being stable in respect to ripening. Thus finite and 
ripened islands coexist in the R2 phase. 
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VW Phase: For large misfit (e > £2) and for small cov- 
erages (H < H C4 (e)), the free energy has its minima at 
ri2 = H and n\ — 0, indicating that all the deposited 
material is accumulated in finite islands. Due to the 
large misfit, in this phase the wetting film is absent and 
the islands form directly on the substrate, similar to the 
Volmer- Weber (VW) growth mode. 

SK2 phase: For €2 < e < £3, increasing H, at H C4 (e), 
we reach the SK2 phase. The behavior of the system is 
different from the SK^ growth mode: at the H CI bound- 
ary we have islands formed in the VW mode. As Fig. 
2b indicates, in the SK2 phase the island density and the 
island size remain unchanged, and a wetting film starts 
forming. This process continues until a full monolayer 
is completed, at which point we enter the SKi phase. 
In contrast with the SKi phase, in the SK2 phase the 
formation of further islands is suppressed until the one 
monolayer thick wetting layer is completed. 

i? 3 Phase: In this phase, present for e > €3 and for 
H > H CA1 the free energy has its minima at n\ — and 
< n 2 < H , indicating the formation of ripened islands. 
The formation of stable islands is suppressed, and all 
the material deposited after H C3 contributes only to the 
ripened islands, coexisting with the stable islands formed 
in the VW growth mode. However, in contrast to R2, in 
R3 a wetting film is absent. 

Comparison with experiments — A quantitative com- 
parison of the phase diagram with the experiments re- 
quires the knowledge of the material constants in (|l])-(||), 
which determine the precise value of ei, £2, £3, and the 
location of the lines in the phase diagram. However, the 
topology of the phase diagram is material independent, as 
long as the SK phase is supported by the system. This ro- 
bustness of the phase diagram implies that in equilibrium 
these are the only phases supported by the free energy 
(@H§- 

During ripening, when the island size reaches a certain 
critical size, dislocation formation relaxes the strain en- 
ergy of the islands, allowing the fast growth of dislocated 
islands. Furthermore, many experiments were done at 
large enough flux so that one suspects that equilibrium 
has not been reached yet. Finally, nucleation barriers 
might slow down the convergence to an equilibrium state, 
trapping the system in metastable states || . 

Our analysis indicates that for strain induced island 
formation there is a critical strain, e±, such that for any 
e > ei stable islands are possible. A second important 
consequence is that for any e, for large enough coverage 
ripening will occur. Thus in order to obtain stable is- 
lands, H must not exceed the boundary of the ripening 
phase. This result can provide an efficient test of our pre- 
dictions: in systems where ripening has been observed, 
we predict that ripening can be avoided by choosing a 
smaller coverage H . 

The formation of the pseudomorphic wetting layer for 
small H and e has been documented in various sys- 
tems, being a general feature of strained layer forma- 
tion. Detailed measurements on InAs/GaAs have shown 



that the transition from the FM to the SKi phase occurs 
at H Cl ~1.7ML, independent of the deposition rate [||, 
indicating that its origin is thermodynamic rather than 
dynamic |icj] . Furthermore, recent investigations have 
measured the strain dependence of H Cl : results on GeSi 
grown on Si and AlInAs grown on AlGaAs have indicated 
that the critical wetting layer thickness decreases with 
increasing misfit pjf| , in agreement with the decreasing 
tendency of the H Cl phase boundary (Fig. 1). 

After the critical thickness has been reached, rapid for- 
mation of uniform islands is observed Q . Studies of InAs 
grown on GaAs indicate that near H Cl the island density 
increases as p ~ (H — H c (e)) J p2[ , signaling a second or- 
der phase transition in the system. Furthermore, we find 
that in the close vicinity of H Cl we have 7=1, i.e. the is- 
land density increases linearly with [H — H Cl ). However, 
for large (H — H Cl ) the island- island interaction leads to 
sub-linear increase in the density. Indeed, Miller et al. 
|]l2| found that after stopping deposition the system had 
a transient regime, after which it equilibrated. The equili- 
brated island density increased linearly with the coverage, 
in agreement with our prediction 7 = 1. Furthermore, a 
linear expression provides an excellent fit near H Cl to the 
data of Ref. Q as well. 

We find that unlike p, the equilibrium island size does 
not increase continuously near H Cl , but it jumps discon- 
tinuously from zero to xo(e, H Cl ). This is again in agree- 
ment with the experiments, since once islands form, they 
reach a well defined size and small islands are rather rare 
El|3[. Also, the experiments indicate that while an in- 
creasing H does modify the equilibrium island size, this 
change is not significant, but most of the newly deposited 
material contributes to the formation of new islands [|IJ, 
again in agreement with slowly changing xq and rapidly 
increasing p in Fig. 2. 

Finally, the phase diagram indicates that the stability 
of the islands depends on the coverage: independent of 
e, for large H ripening should take place in the system. 
Indeed, for CdSe grown on ZnMnSe repeated AFM scans 
of the same sample made at 48 hour intervals indicate 
that some islands ripen at the expense of others, and the 
overall island density decreases with time [jl4|. However, 
the fact that the stable islands do not coexist with the 
ripened ones, as expected in the R2 and R3 phases, sug- 
gests that these experiments were performed either at 
the FM and Ri phase boundary, or dislocations relax the 
ripened islands. The coexistence of stable and ripened 
islands is documented for Ge grown on Si jig ], that, to- 
gether with the evidence of a wetting film in this system, 
indicate that these experiments are done at the border 
of the SKi and R2 phases, and consequently stable SK 
islands are allowed for smaller coverages, assuming that 
dislocations are not the sole origin of the ripened islands. 

We have benefited from useful comments and discus- 
sions with J. K. Furdyna and J. Tersoff. 



3 



[1] For a recent review see W. Seifert, et al, J. Progr. Crys- 
tal Growth Charact. Mater. 33 423 (1996); P.M. Petroff, 
and G. Medeiros-Ribeiro, MRS Bulletin 21 (No 4), 50 

(1996) . 

[2] J.-M. Gerard, in Confined Electrons and Photons, Ed. 
E. Burstein and C. Weisbuch (Plenum Press, New York, 
1995). 

[3] J. Tersoff, Phys. Rev. B 43, 9377 (1991); C. Roland and 

G. H. Gilmer, Phys. Rev. B 47, 16286 (1993). 
[4] A.-L. Barabasi, Appl. Phys. Lett. 70, 2565 (1997); B. G. 

Orr et al, Europhys. Lett. 19, 33 (1992). 
[5] V. A. Shchukin, et al, Phys. Rev. Lett. 75, 2968 (1995). 
[6] L. D. Landau and E. M. Lifshitz, Theory of Elasticity 

(Pergamon Press, Oxford, 1986). 
[7] V. I. Marchenko, Sov. Phys. JETP 54(3), 605 (1981). 
[8] J. M. Rickman and D. J. Srolovitz, Surf. Sci., 284, 211 

(1993). 

[9] D. Jesson et al, Phys. Rev. Lett. 77, 1330 (1996). 
[10] I. Daruka and A.-L. Barabasi, Phys. Rev. Lett. 78, 3027 

(1997) . 

[11] G. Abstreiter, et al, Semic. Sci. Tech. 11, 1521 (1996). 
[12] D. Leonard, et al, Phys. Rev. B 50, 11687 (1994); N. P. 

Kobayashi, et al, Appl. Phys. Lett. 68, 3299 (1996); M. 

S. Miller, et. al, Sol. St. Elec. 40, 609 (1996). 
[13] T. I. Kamins, E. C. Carr, R. S. Williams, and S. J. Ros- 

ner, J. Appl. Phys. 81, 211 (1997). 



[14] S.H. Xin, et al, Appl. Phys. Lett. 69, 3884 (1997). 
[15] R. Leon, et al, Appl. Phys. Lett. 67, 521 (1995). 



FIG. 1. Equilibrium phase diagram in function of the cov- 
erage H and misfit e. The small panels on the top and the 
bottom illustrate the morphology of the surface in the six 
growth modes. The small empty islands indicate the presence 
of stable islands, while the large shaded one refer to ripened 
islands. The phases are separated by the following phase 
boundary lines: H ci (e): FM-Ri, FM-SKi; H C2 (e): SKi-R 2 ; 
H C3 {e): SK2-SK1; # C4 (e): VW-SK 2 , VW-R3. The parame- 
ters used to obtain the phase diagram are a = 1, C — 40i5o, 
$aa = E , $as = 1.27£ , 9 = 0.7, p = 4.9, 7 = 0.3, b = 10. 

FIG. 2. (a) Wetting film thickness (rii), island coverage 
(112) (top), island size (xo), and island density (p) (bottom) 
as a function of H for e = 0.08. At H C1 there is a transition 
from the FM to the SKi phase, the island size jumping dis- 
continuously. In the R2 phase, present for H > H C2 , ripening 
takes place; (b) Same as (a) but for e = 0.12. At H C4 there is 
a transition from the VW to the SK2 phase followed by the 
SKi phase at H C3 . And finally, at H C2 the system reaches the 
R2 phase. 
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